New collective mode due to collisional coupling 
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Starting from a nonmarkovian conserving relaxation time approximation for collisions we derive 
coupled dispersion relations for asymmetric nuclear matter. The isovector and isoscalar modes 
are coupled due to asymmetric nuclear meanfield acting on neutrons and protons differently. A 
further coupling is observed by collisional correlations. The latter one leads to the appearance of 
a new soft mode besides isoscalar and isovector modes in the system. We suggest that this mode 
might be observable in asymmetric systems. This soft mode approaches the isovector mode for 
high temperatures. At the same time the isovector mode remains finite and approaches a constant 
value at higher temperatures showing a transition from zero sound like damping to first sound. The 
damping of the new soft mode is first sound like at all temperatures. 
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The investigation of collective excitations in asym- 
metric nuclear matter is of current interest for experi- 
ments with nuclei far from (3- stability, jl| and citations 
therein. We consider here a Fermi gas model consisting 
of a number of different species (neutrons, protons, etc) 
interacting with the own specie and with other ones. The 
interaction between different sorts of particles is impor- 
tant to consider if we want to include friction between 
different streams of isospin components. Especially the 
isospin current may not be conserved by this way. We 
neglect explicitly shell effects and concentrate only on 
bulk matter properties. Let us start with a set of cou- 
pled quantum kinetic equations for the reduced density 
operator p a for the specie a 



dtPa(t) = i[p a ,£a + U(t) a ] - 
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dt 



,p a {t')-p b {t') 

T ab (t-t>) 



(1) 



where £ = V 2 /2m denotes the kinetic energy operator 
and U the mean field operator and the external field 
which is assumed to be a nonlinear function of the den- 
sity. We have approximated the collision integral by a 
non-Markovian relaxation time The memory effects 
condensed in the frequency dependent relaxation time 
turned out to be necessary to reproduce damping of zero 
sound (^||. It accounts for the fact that during a two 
particle collision a collective mode can couple to the scat- 
tering process. Consequently, the dynamical relaxation 
time represents the physical content of a hidden three 
particle process and is equivalent to the memory effects. 

We have further assumed the relaxation with respect to 
the local equilibrium p b of any specie in the system. The 
relaxation of the actual distribution of specie a is driven 
by the local equilibrium of all the other components. The 
cross coupling is realized by nondiagonal relaxation times 
Tab- We specify the local equilibrium by a small deviation 



of the chemical potential of specie a 
the global equilibrium 



compared with 



< k\p a \k' >= f a (k)S kk , - fa{k J fa ZU va(k ~ k') (2) 

e a(. fc ) € a\' i ) 

with < k\£ \k' >= e(fc). The equilibrium distributions 

< k\p° a \k' >— f a (k)Skk' are the corresponding Fermi 
functions with chemical potential p a and the normaliza- 
tion to density n a — 2 f j0fjsfa(p)- The local equilibrium 
specified by Sp a is determined if we impose the density 
balance to be fulfilled separately for each specie current 
J a which reads in Wigner coordinates k = p + q/2 and 
k'=p-q/2 



u)5n a (q,uj) = qdJ a (q,uj). 



(3) 



From this equation we derive the following matrix equa- 
tion for Sp a 



n Q (g,0)^ a = ^^| ^ {B b 5(q)+ 5 -^ 



(4) 



where B a = £ 

b 
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and the partial polarization function of specie a is 



!!„(<?, w) = 2 



dp /a(p+§)-/a(p-f) 

(2tt)3 e a (p+l)-e a {p~±)-Lo- 



(5) 



(6) 



The factor 2 in front of the integral accounts for the spin 
degeneracy. Eq. (^) generalizes eq. 6 of 0] to multicom- 
ponent systems. 

Now we linearize the kinetic equation (Q) around the 
total equilibrium p° a which we express as 



1 



Spa = Pa - Pa = Pa - Pb + Pb 



Pa- 



iJ) 



This provides us with the density variation 5n due to 
an external field U cxt which is connected with the po- 
larization of the system via H5n = U cxt . The obtained 
polarization function has now a matrix structure. The 
poles of the polarization function represent the collective 
modes in the system. For a two component system, e.g. 
neutrons with density n n and protons with density n p , 
we have the density variation of the mean field 



SU a = a an Sn n + a ap 5n p 



(8) 



and we obtain for the collective modes the dispersion 
relation 



(l-n^OnnXl-II^app) 

-{D np + n„ a np )(D pn + n a pn ) = 0. 



(9) 



The generalization of the Mermin polarization function 
Q to a multicomponent system is derived with the inclu- 
sion of nonmarkovian (frequency dependent) relaxation 
times as 
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n (w + i) 



i - 



U)T a -\-i 



(1 - C aa ) 



(10) 



An additional coupling in (^]) occurs due to asymmetry 
and collisions 



D r , 



(ii) 



The D pn are given by interchanging sort indices. The 
matrix C ao is expressed as 



H c ((u; + ^) 

n c (o) 



(12) 
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The term D np is vanishing for symmetric matter, i.e. for 
equal densities of species as well as for vanishing collision 
integral r — > oo. Therefore we call this term asymmetry 
coupling term further on. 

The dispersion relation (|^) is similar to the one derived 
recently in |lj if we neglect the collisional coupling D np . 
The latter one has also been discussed in plasma two- 
stream instabilities ||. Here we present a more general 
dispersion including correlational coupling. The polar- 
ization function includes collisions within a conserving 
approximation ||. 

Before we apply this dispersion relation to nuclear mat- 
ter let us consider a special case. We assume symmetric 



a pp — cti and a 



nuclear mean fields a n 
and neglecting the collisions we obtain 



np 



x-pn 



1 — (ai ± aa) 



n( w + i) 



-M 1 

U)T-\-l v 



n(o) > 



0. 



a 2 



(13) 



These are the decoupled dispersion relations for the 
isovector mode a% — a 2 and the isoscalar mode a\ + a 2 . 

Summarizing, we see that there appear two kinds of 
coupling: (i) The coupling of modes between isovector 
and isoscalar ones due to different mean fields and (ii) 
an explicit correlational coupling of asymmetric nuclear 
matter due to collisional correlations, which is condensed 
in D np . We have presented a general dispersion relation 
for the multicomponent system including known special 
cases. 

In the following we will apply this expression for the 
damping of giant dipole resonances in asymmetric nu- 
clear matter. We assume a wave vector for giant dipole 
resonances corresponding to the Steinwedel and Jensen 
JtJ model q — ^ with the nuclear radius R = 1.13 A 1 / 3 
fm connected to mass number A. This wave vector has 
very low values compared with the Fermi wave vector. 
Therefore it allows us to expand the Mermin polariza- 
tion function with respect to small qv c /oJ ratios where v c 
is the sound velocity. The frequency dependence of the 
dynamical (memory) relaxation times is derived using a 
Sommerfeld expansion S 
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Tafc(w) T afc (0) 
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4 VttT/ 



(14) 



for a, b neutrons or protons respectively. The markovian 
relaxation time was given in terms of the cross section 
a a b between specie a and b as r^ 1 = ^-a a i,T 2 . 

The dispersion relation (|j|) takes then the form of a 
polynomial of tenth (six) order corresponding to the in- 
clusion of memory (in)dependent relaxation times via 



= ( w ( w + — ) - cl n q 2 )(uj(u + —)- c 2 pp q 2 ) 




with the definition (||) and (14). Here the partial sound 
velocities are 



2 n a {na) 

c ab = a ab 

m 

-2 1 E„ — E„ 
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and denoting the Fermi energy with 

,/3/ 2 (e*?) 



E a =T- 



/i/ 2 (e^) 



(16) 
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We observe that the new collisional coupling represented 
by c 2 b vanishes if cither the diagonal friction t^ 1 or the 



2 



nondiagonal friction t^ 1 vanishes or the system is sym- 
metric E n = E p . This underlines that we have a new 
coupling due to collisional correlations and asymmetry. 

We use as an illustrative example the following mean 
field parameterization of Vautherin P,|lC|] 



+j{(n n +n p ) 2 -n 2 a ) 



(18) 



with n a = n n ,n p the density of neutrons or protons, re- 
spectively. The corresponding mean field deviations can 
easily be computed via a a t, — dU a /dni,. The Coulomb 
interaction leads to an additional contribution for the 
proton meanficld 



(19) 



The here used model parameters reproduce the 
Weizsacker formula 



E a 2 a 3 Z 2 2 



(20) 



In figure [l] we plot the occurring modes for a weak 
asymmetric case of 11 B. We observe that the isovector 
mode is decreasing with increasing excitation energy and 
vanishes at about 5 MeV. The corresponding isoscalar 
mode is already vanishing at 4 MeV. In general all en- 
ergies appear as symmetric solutions with positive and 
negative energy where negative energy solutions are ruled 
out as unphysical. The corresponding damping is de- 
generated up to the point of vanishing energy. Above 
this temperature the damping of isoscalar and isovector 
modes become twofolded. 

11 B ( 5=0.1 ) 
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by the volume energy oi = 15.68 MeV, Coulomb energy 
a 3 = 0.717 MeV and the symmetry energy 04 = 28.1 
MeV with the asymmetry parameter 



(21) 



In figure y\ we have plotted the solution of the dispersion 
relation (P) for 11 Be versus excitation energy. The ki- 
netic energy is linked to a temperature within the Fermi 
liquid model via Sommerfeld expansion 




E e „ [ MeV ] 



E 3 /(l + <5) 5/3 + (l-5) 5 / 3 



+ T2^T f 



T , 2 (1 + ^) 1/3 + (1-<5) 1/3 



(22) 



This connection between temperature and excitation en- 
ergy is only valid for a continuous Fermi liquid model. 
For the small nuclei like n Be, the concept of temper- 
ature is questionable. Some improvement one can ob- 
tain by the definition of temperature via the logarithmic 
derivative of the density of states JllJ 



T~ 



1 dp 



pdE c 



A 



4e f E e 



a/2 



(23) 



which provides i? cx w j(E/A) in comparison with ( |2^ ) 
for small temperatures. We use this temperature to 
demonstrate possible collective bulk features in an ex- 
ploratory sense. Of course, the surface energy and shell 
effects cannot be neglected for realistic calculations. 



FIG. 1. The centroid energy (above) and the damping 
width (below) of collective modes vs. excitation energy 
Ecx = \{E/A) (above Fermi energy) found for 11 B situation. 
The isovector mode (rhombus) vanishes at about 5 MeV while 
the corresponding isoscalar mode (triangles) already disap- 
pears at 4 MeV. Connected with the disappearing the corre- 
sponding damping becomes twofolded. This is explained by 
the symmetric negative modes plotted as thin lines which dis- 
appear as well. A soft mode (squares) arises and is connected 
with an continuous increasing damping width. The inset gives 
the instable modes separately. 

Besides the standard isovector and isoscalar modes we 
observe a build up of a very soft mode with a centroid 
energy around 1 MeV. This mode appears due to the 
collisional coupling c 2 ab of (|l6|). When we turn off the 
relaxation times, i.e. the collision integral, this mode is 
vanishing as well as in symmetric nuclear matter, see dis- 
cussion after ([!(]). It shows that this mode appears due 
to collisional coupling of isovector and isoscalar modes. 
The corresponding damping of the crossed mode is con- 



3 



tinuously increasing with temperature. 

Due to the consideration of memory (frequency depen- 
dent) relaxation times we observe further an unstable 
mode shown in the inset of figure |l|. If we neglect this 
memory effect we would not observe this instable modes 
but the soft mode described above remains. 
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FIG. 2. The centroid energy (above) and the damping 
width (below) of collective modes vs. excitation energy (above 
Fermi energy) found for 11 Be situation with the same set up as 
figure [l]. The isovector mode (rhombus) decreases and remain 
constant at about 5 MeV while the corresponding isoscalar 
mode (triangles) already disappears at 4 MeV. A soft mode 
(squares) appears with higher centroid energy than in 11 B of 
figure [l] while the damping width is the same. 

Looking at a more asymmetric situation of 11 Be in fig- 
ure |^ we see that the centroid energy of the third mode 
becomes larger while the damping remains the same. A 
new feature appears for the isovector mode. While in 
weak asymmetric systems the mode disappear sharply 
at a certain temperature, we see now a decrease and a 
convergence towards the new soft mode. This suggest 
that the new soft mode should be of isovector character. 
This behavior of the isovector mode is connected with no 
transition to a twofolded damping we have seen for the 
more symmetric case. Instead the damping of isovector 
modes show a transition from a T 2 temperature behavior 
typical for zero sound damping towards a E 1 ^ 4 behavior 
typical for first sound damping. We like to pronounce 
that in symmetric nuclei we did not observe this tran- 
sition because the mode itself, i.e. the centroid energy 
is vanishing at the point where this transition behavior 



of the damping occurs. In asymmetric nuclear matter it 
should be possible to observe this transition from zero 
to first sound because the energy of the isovector mode 
remains finite. The twofold damping of isoscalar mode 
is modified with respect to the symmetric case but not 
removed. 

One may argue whether this third mode can really ap- 
pear in the system. A simple consideration may con- 
vince us about the possible existence of such mode. Let 
us assume a coupled set of two type of harmonic oscilla- 
tors (neutrons and protons) interacting between the same 
sort of particles with strength k n and k p , respectively 
and between different sorts with k np . Let us choose for 
simplicity only two neutrons and two protons. Then we 
obtain the coupled system of harmonic oscillators with 
frequencies ui„ = k n /m, uj p — k p /m and w 2 ^ = k np /m. 
The solution yields three basic modes in the system, 
i.e. w 2 = 2(w 2 + w 2 p ), 4w 2 p , 2{uj 2 p + uj 2 lp ). If we neglect 
the different coupling between neutrons and protons io np 
we only obtain two modes analogously to isovector and 
isoscalar ones. We see that the coupling between neu- 
trons and protons can lead to the appearance of a third 
mode. 



"Be, B, C, N ( 5=0.3, 0.1, -0.1, -0.3 ) 




2 4 6 8 10 

E m I MeV ] 




E e „ [ MeV ] 



FIG. 3. The centroid energy (above) and the damping 
width (below) of collective modes vs. excitation energy (above 
Fermi energy) for the new soft mode for different isobaric 
states of 11 Be. The symmetric nuclei lx Be (circles) and 11 N 
(rhombus) are almost equal and three times larger in the cen- 
troid energy than the more symmetric cases 11 B (triangles) 
and 11 C (squares). 

We conclude that due to collisional coupling there 
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should be a third soft collective mode observable besides 
the isovector and isoscalar modes. To study further these 
modes may open an exciting access to correlation effects 
in nuclear matter. In figure |^ we give the predictions for 
the isobar states of u Be. We find that for U B and 11 C 
the modes are by a factor three lower. Interesting to re- 
mark that the damping width is found to be independent 
of the asymmetry. 

Let us now compare the found new mode with the 
experimental evidence. There are some hints for a soft 
mode in 11 Be Hj. The authors have observed a low 
lying structure at around 6 MeV excitation energy with a 
damping of around 1 MeV which has not been reproduced 
yet even within refined coupled channel calculations [fl3[ . 
A standard explanation would give as the origin a weekly 
bounded single particle neutron orbital. The observed 
broad structure at 6 MeV might be explained possibly as 
the here presented coupled mode. The centroid energy 
as well as damping width at least seem to suggest this 
interpretation. 

Again, we like to point out that the presented results 
are limited to a pure liquid drop model. We have ne- 
glected shell effects and surface effects which are impor- 
tant for realistic calculations in small systems like 11 Be. 

To summarize we have observed that due to correla- 
tional coupling there can exist a new mode which appears 
besides isovector and isoscalar modes in asymmetric nu- 
clear matter due to collisions. We suggest that this mode 
may be possible to observe as a soft collective excitation 
in asymmetric systems. The transition from zero sound 
damping to first sound damping behavior should become 
possible to observe for isovector modes since they do not 
vanish at this transition temperature like in symmetric 
matter. 
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